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ABSTRACT 

The analysis of eight EXOSAT X-ray hghtcurves of six active galactic nuclei by non- 
linear prediction methods indicates that the observed variability is truly stochastic 
and is not caused by deterministic chaos. This result favors X-ray emission models 
with mutliple centers of production of hot, possibly relativistic electrons. 
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1 INTRODUCTION 

X-ray emission is clearly one of the key points in under- 
stand the nature of radio quiet AGNs. It comprises a sub- 
stantial fraction of the total bolometric luminosity and is 
strongly variable on time scales as short as hours (see, e.g. 
Mushotzky, Done & Pounds 1993). Any AGN model must 
therefore explain both the level of X-ray emission and its 
rapid variability. 

The longest (up to 2 x 10^ s) uniterrupted light curves 
for short time scale variability studies in the 1-8 keV range 
have been provided by the EXOSAT due to its unusual orbit. 
Longer time scale variability has been studied by combining 
data from different satellites over several years. 

Several characteristics of X-ray variability of AGN have 
emerged in these studies. The X-ray flux is strongly vari- 
able (e.g. McHardy 1988, Green, McHardy & Lehto 1993) 
in most sources. Power spectra are featureless and have a 
power-law shape S o; , with q ~ — 1 over a broad range of 
frequencies (~ 10"^ - 10"^ Hz) (McHardy & Czerny 1987, 
Lawrence et al. 1987, Lehto et al. 1991; for a review see, e.g. 
McHardy 1990). The only features seen in these spectra 
are the expected flattening of the spectrum at timescales of 
years, which is claimed for NGC 5506, and the signs of quasi- 
periodicity in NGC 5548 (Papadakis & Lawrence, 1993) and 
NGC 4051 (Papadakis & Lawrence, 1995). Note that the 
claimed periodicity in NGC 6814 (Mittaz & Branduardi- 
Raymont 1989; Done et al. 1992) was subsequently discov- 
ered to be the property of the galactic X-ray source con- 
taminating the X-ray measurements in the direction of this 
AGN (Madejski et al. 1993). 

As the primary emission strongly dominates in the EX- 
OSAT energy range we can assume that the conclusions ap- 
ply directly to the primary power-law component of the X- 
ray emission. 

An observed featureless power spectrum is expected in 
two qualitatively different dynamical systems: 



(A) a system consists of a large number of uncorre- 
lated (weakly correlated) elements varying (almost) inde- 
pendently 

(B) global evolution of the system is described by three 
or more non-linear differential equations which shows deter- 
ministic chaos behaviour. 

These two interpretations lead to two different physical 
pictures. 

Case (A) can be interpreted, for example, as randomly 
generated shocks in outflowing (or accreting) gas, or ran- 
domly appearing reconnections of large scale magnetic field 
lines, in close analogy to phenomena in solar corona. In this 
case the emitting medium is clearly clumpy and the emission 
itself can be modelled as a shot noise (e.g. Lehto 1989). 

Case (B) could correspond to coherent variability of a 
single cloud of electrons, e.g. some generalization of opti- 
cally thin geometrically thick disk, spherically accreting hot 
(possibly pair dominated) plasma or instability of a single 
shock at the base of jet forming region. 

Sugihara & May (1990) suggessted that nonlinear pre- 
diction methods are (in principle) capable of distinguishing 
between cases (A) or (B). 

The methods (Tsonis & Eisner 1992) are based on the 
study of the decay of short term correlations in the light 
curve. Exponential decay of the correlation is character- 
istic for a chaotic signal with a strange attractor while a 
power law decay is expected from shot noise (more precisely, 
from the so called fractional brownian motion representing 
stochastic phenomena with power law power spectra) . 

We apply this new method to a number of long EX- 
OSAT observations of Seyfert galaxies to which the correla- 
tion integral method has been applied previously by Lehto, 
Czerny & McHardy (1993; hereafter LCMH). 



2 B. Czerny & H.J. Lehto 



2 METHOD 

Wc apply the non-linear prediction methods of Tsonis and 
Eisner (1992) to investigate the characteristics of X-ray vari- 
ability of several AGNs that show boardband l/.f" type 
power spectra. The essence in these methods is to break up 
the data into short strings of consecutive data points and 
then to analyse the statistical properties of these strings. 
They are far more efficient than the Fourier methods in 
studying the behaviour of a system with short term cor- 
relations. 

The method compares shapes of two data strings of 
length d. Mathematical expression of that procedure is the 
calculation of the distance between the two strings, or vec- 
tors X = {xi,X2, .., Xd) and Y = (j/i, j/2, ■■■^Ud)- Usually the 
distance r, or the norm, is defined as the L2 norm of the 
difference vector: 



in close analogy to St of two curves. We adopt this 
definition in our paper. Other norms can also be used, for 
example, the Chebychev or Loo norm, 

r* = max (jxi — (2) 

i=l,d 

The aim of the nonlinear prediction method is used in 
an attempt to predict the position of the point separated by 
time interval T from the sequence Y on the basis of known 
library data set Xk- 

This consists of three steps: (i) identiiying the se- 
quences from the library {Xk} most similar to the one be- 
ing considered (Y) (by minimizing the norm between) (ii) 
identifying values the strings in the library would get af- 
ter a timestep T (iii) calculating the predicted value using 
the values from (ii) (iv) and finally correlating it with the 
observed value. 

We have used three different methods for step (iii): 

• Nearest Point Method. Identify the nearest X as the best 
predictor sequence. Use that sequence alone to predict Y 
after at time interval T. 

• Linear Interpolation Method. Identify the d + 1 closest 
sequences. Use linear interpolation in d dimensional space 
to obtain a predicted value after a time interval T from the 
predictor sequences. 

• Exponential Weight Method. Identify d -I- 1 closest se- 
quences. Predict Y at time interval T by using the observed 
values of the X sequences, but now interpolating them with 
a weight of exp(-r) for each predictor value, where r is the 
distance between the given sequence and Y. 

We divide the data sequence into two halves. The first 
part forms the library and the second one is used to test the 
quality of prediction. Since we actually know the values, we 
have a set of pairs {pk,fk) consisting of the predicted and 
known value for step (iv). 

The quality of prediction can be evaluated by calculat- 
ing the correlation coefficient R between the series of pre- 
dicted pk and known fk values: 

R{T) = {(pkfk)-(j:>k){fk)}/T^ (3) 



where {•) denote the average of the quanity over the series 
and a, the standard deviation of the mean value of the time 
series. 

If the data is dominated by short term correlations the 
quality of prediction decreases with time step T. Tsonis & 
Eisner (1992) shows that for, for chaotic signal (in a sense 
of deterministic chaos) 

1 - R{T) DC exp {2KT) (4) 
and in the case of fractal brownian motion (shot noise) 

1 - R{T) oc T^" (5) 

Here K is the Kolmogorov entropy (Kolmogorov 1959) 
which is the sum of the positive Lyapunov exponents (see 
e.g. Schuster 1989) and H is related to the slope a of the 
power density spectrum through the relation a = 2H + 1. 

By correlating log(l — R) againt T and against log(T) 
we obtain Ra and Rs, two further correlation coreffincients 
which tell us how well each of the two predictive modelst 
decribes the data. We identify the preferred system by se- 
lecting the system with the larger of Rs and Rd. 

3 RESULTS 

3.1 Tests of the method using NGC 4051 ME 

We tested technical details such as the choice of length 
string, d, the binning, the method of prediction and the time 
range of prediction on the EXOSAT ME X-ray lightcurve of 
NGC 4051. 

Binning the data is the standard approach to increase 
the signal to noise ratio. This is done at an expense of the 
available dynamical range. 

We used generally 200 second bins, corresponding to 
S/N ratio of 6 in NGC 4051. We made also some further 
tests with 50 second bins. 

The value of the correlation coefficient depends on the 
assumed length of the string d. Therefore we first calculated 
this coefficient for d ranging from 1 to 20 and the prediction 
time step T = 1 (one bin). Wc show this dependence in fig.l 
using the three different prediction methods. The linear 
interpolation method was used only up to d = 4, as the 
search for the smallest simplex containing a given vector 
(c.f. Sugihara and May 1990) is time consuming. 

We can see that the third method, i.e. the exponential 
weight gives the highest values of the correlation coefficient 
and relatively small wiggles around the systematic trend. 
Therefore, in all the following computations, we use the ex- 
ponential weight method to compute the predicted values. 

The value dmax corresponding to the best prediction is 
picked up for the purpose of further calculations for larger 
T. 

We also tested the prediction for the same time interval 
(i.e. 200 s) but with bins of 50 seconds. The highest values 
for the correlation coefficient were met for strings approxi- 
mately longer by a factor 4, as expected, but the values of 
the correlation coefficient was lower. It was also lower for 
the prediction for a single time step (i.e. 50 s). Too strong 
binning will lead to a loss of information about the corre- 
lations in the system. Assuming time step 800 s (binning 
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Figure 1. The dependence of the correlation coefBcient on di- 
mension of the string used for prediction. NGC 4051 ME Ught 
curve is used as an example. Prediction time step equal 1 bin (200 
s). Three prediction methods were used; nearest point (squares), 
linear interpolation (stars) and exponential weight (triangles). 

again by a factor 4) results also in the decrease of the corre- 
lation coefficient of the prediction for a single time step (i.e. 
800 s). Our initial 200 second binning appears therefore 
reasonable. 

Next wc calculate the prediction quality as a function 
of prediction time. We use 200 second bins and strings of 9 
points. Figure 2 shows the result of plotting log(l — R) vs. 
log(r) and log(l - R) vs. T 

The distinction between the two cases is not very pro- 
found, but a clear tendency is seen to favor the log-log de- 
pendence, which corresponds to the fractal brownian type 
shot noise variability. The correlation coefficient Rs is 0.998, 
and the coefficient Rd is equal 0.958. 

We test the dependence of the measured correlation de- 
cay on the choice of the string length. In Fig. 3 wc show 
the same plot as in Fig 2b but derived using different string 
length dmax ■ We see that for the string length dmax = 9 the 
plot is very smooth whilst shorter length values introduce 
more wiggling. However, in all cases the correlation coef- 
ficient Rs is greater than Rd and the slopes of the curves 
are similar. Therefore the use of the string length giving 
the best prediction for a single time step seems to be well 
justified. 

We may be tempted to obtain more significant result by 
increasing the dynamical range of the study, that is making 
predictions for even longer time steps. However, the decay of 
correlations takes us into the range dominated by the data 
noise after some 15 steps, i.e. interval of 3000 seconds for 
this source. 

Since the shot noise is favored the slope, 2H, of the 
curve log(l-R) vs. log(T) should be related to the slope of 
the power spectrum, o, through the relation a = 2H + 1. 
As 2H — 0.42 and the power spectrum slope determined 
by Papadakis & Lawrence (1995) for this lightcurve is a = 
1.461q Jg, the two numbers are in agreement. 
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Figure 2. The quality of prediction (see eq. 4 and 5) as a function 
prediction time step for exponential method. The decay of corre- 
lation, i.e. decrease of R is clearly visible. The rate of this decay 
is better represented by a straight line on plot log(l-R) vs. log(T) 
(lower panel) than log(l-R) vs. T (upper panel). Best power-law 
fit to the decay law is shown with continuous line whilst best ex- 
ponential fit to decay is shown with dashed line on both panels. 
The respective correlation coefficients for these fits Rs and Rd 
are equal to 0.998 and 0.958, thus stochastic interpretation of the 
lightcurves is clearly favored. 



3.2 Simulated lightcurves 

In order to test the sensitivity of the method we use two 
computer generated light curves of the equal lengths. The 
first one is the Lorenz attractor (Lorenz 1963), an example 
of a source exhibiting deterministic chaos. We use a time 
step 0.02 and the coordinate X to represent the light curve 
as in our previous paper (LCMH). The correlation decay 
plot for this lightcurve is shown in Fig. 4. 
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Figure 3. The dependence of the quahty of prediction on the 
adopted length of the string. The values used are: 1 (filled train- 
gles), 4 (open traingles), 9 (filled squares) and 10 (open squares). 
All values favour stochastic interpretation of the lightcurve. The 
value 9, giving the highest correlation in prediction for a single 
step is adopted in further computations for this lightcurve. 
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Figure 4. The decay of correlations as a function of prediction 
time step for a simulated light curve formed from the X-coorinate 
of the Lorenz attractor. Exponential decay (dashed line) is fa- 
vored (see eq. 4) over the power law trend (continuous line). 
This is expected because this case is an example of a simulated 
light curve caused by deterministic chaos. 

As expected, the character of variability is better de- 
scribed by exponential decay of correlations, and the coef- 
ficient Rs is smaller than Rd (0.949 and 0.986 respectively, 
for the length of lightcurve 2000 points) . 

The second computer generated lightcurve is an exam- 
ple of shot noise (Lehto 1989). Again, the nonlinear predic- 
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Figure 5. The decay of correlations as a function of predic- 
tion time step for computer generated example of stochastic light 
curve; as expected, this time power law trend (continuous line) is 
favored (see eq. 5) over the exponential law (dashed line). 

tion method unveils correctly the nature of the light curve 
{R, = 0.99, Rd = 0.92). 

3.3 EXOSAT background 

The effects of background subtraction may, in principle, in- 
fluence the results of a variability study (see e.g. Fiore & 
Massaro 1989). Fig. 6a shows the decay of correlations for 
NGC 4051 ME lightcurve (background subtracted) and the 
background itself. Clearly the behaviour of the two is dif- 
ferent. The level of correlations in the background is quite 
high and it almost does not decay with the increase of the 
prediction time. The formal fit slightly favors the interpre- 
tation of deterministic chaos, as the Rd and Rs coefficients 
for the ME background during the observation of NGC 4051 
are equal to 0.904 and 0.878, respectively. Such long range 
correlations may reflect some systematic trends of unknown 
nature. However, it is encouraging, that correlations present 
in the signal are markedly different from those present in the 
background. 

In the case of LE EXOSAT lightcurve of NGC 4051 
(Fig. 6b), the decay of correlations for the source is smooth 
and nicely follows a straight line on log-log diagram; the 
behaviour of the background is badly approximated by any 
of the two laws. It shows again that the correlations are 
not just artifacts of the background subtraction but the real 
property of the source. 

3.4 EXOSAT lightcurves 

The nonlinear prediction method was applied to a number 
of long EXOSAT X-ray lightcurves. The data was binned 
in 200 s. The exponential weight method was used for pre- 
diction. We excluded the light curve of NGC 6814 because 
it is strongly contaminated by the X-ray emission from a 
variable galactic source (Madejski et al. 1993). 
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Table 1. The correlation coefficients for power law and exponen- 
tial decay of correlations in EXOSAT lightcurves. 
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Figure 6. The decay of correlations as a function of predic- 
tion time step for the source liglitcurvc (NGC 4051) (squares) 

and the background (circles) from the two EXOSAT instruments; 
ME (upper panel) and LE (lower panel). Straight lines show the 
best fit to the decay of correlations under the assumption of the 
stochastic nature of the signal. 

Seven light curves arc from EXOSAT ME instrument 
sensitive to 1 - 8 keV radiation. One lightcurve of NGC 
4051 is from EXOSAT LE instrument operating in 0.25 keV 
- 1.5 keV band. 

The correlation coefficients i?s, Rd of the light curves 
are shown in Tabic 1 with some other relevant information. 

The correlation cocfficcicnts of the light curves appear 
to fall into two classes. The same separation is also evident 
if plotting the decay of prediction quality as a function of 
adopted time step. 

Class 1 contains Mkn 335, NGC 4151 and 3C 273. AU 
these sources are characterized by rather low values (\R\ < 
0.2) of both correlation coefficients for time steps greater 



Name 


Obs.time(103 s) 


Rd 




Class 


Mkn335 


67 


0.557 


0.562 


1 


NGC4151(93) 


62 


0.580 


0.753 


1 


NGC4151(95) 


37 


0.544 


0.741 


1 


3C273 


142 


0.557 


0.750 


1 


MCG6-30-15 


178 


0.977 


0.993 


2 


NGC5506 


239 


0.947 


0.984 


2 


NGC4051iB 


202 


0.843 


0.959 


2 


NGC405lMi5 


196 


0.954 


0.999 


2 



than 1 bin although the prediction for one bin is sometimes 
significant - the same effect is observed in the behaviour of 
the background lightcurve (see Sec. 3.2). Visually the four 
light curves appear to show little variation on short times 
scales. On timescales of hours slowly varying trends are 
apparent. 

No conclusion on the character of variability in these 
sources can be drawn on formally fitted coefficients Rs and 
Rd, which are not significant (lower than ~ 0.8). The ratio 
of the two coefficients is in all four light curves is greater 
than 1, which may favor stochastic behaviour. However, the 
same is true for the background lightcurve. It appears that 
the low correlation coefficients can be understood as being 
caused by the dominance of poisson noise to real variability 
of the source on short timescales. 

Class 2 contains NGC 4051 (both ME and LE), MCG 
6-30-15 and NGC 5506. All these sources show a range 
of prediction time steps with high prediction quality and 
the decay of prediction quality is smooth. The correlation 
coefficients Ra and Rd are highly significant, and the ratio 
is larger than 1. 

This result means that (at least in class 2 sources) the 
variability is caused by a truly stochastic phenomenon and 
not by the deterministic chaos. 

The nonlinear prediction method used in his paper is 
basically different from the correlation integral method used 
to these sources in our previous paper (LCMH). Therefore 
it is interesting to note that the classification of the sources 
proposed by us has not changed much. Class 1 sources from 
LCMH (white noise) all found their place in a new class 
1 (sources without good description of short term correla- 
tions). A source classified as class 2b in LCMH (NGC 4151) 
moved to a new class 1. All old class 2 sources (sliot-noisc 
type behaviour) are again in cleiss 2 and the stochastic char- 
acter of their variability found a strong support from the 
nonlinear prediction method. 

The key result is that the lightcurve NGC 4051LB sus- 
pected of chaotic behaviour on the basis of the correlation 
integral method now revealed itself as just another stochas- 
tic phenomenon. 

It is also interesting to note that the character of vari- 
ability of the source NGC 4051 in LE and ME band is very 
similar, the slopes, 2H, for the two lightcurves are equal 
0.43 and 0.42, correspondingly, as if the (at least variable 
fraction of) radiation comes from the same region and the 
same emission mechanism. This conclusion was by no means 
obvious as in many Seyfert galaxies the soft X-ray band is 
dominated by the soft X-ray excess and the presence of an 
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Figure 7. The decay of correlations as a function of the predic- 
tion time step for sources from the Class 1: filled squares - 3C273, 
filled triangles - NGC 4151 in 1983, open squares - NGC 4151 in 
1985, and open triangles - Mkn335. 
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Figure 8. The decay of correlations as a function of the pre- 
diction time step for sources from the Class 2: filled squares - 
NGC 5506, filled triangles - MCG 6-30-15, open squares - NGC 
4051mb, open triangles - NGC 4051_l_b 

excess, albeit not a strong one, was also invoked in the case 
of NGC 4051 (e.g. Walter & Fink 1993). 

What is more, the slope of the power spectrum of LE 
curve calculated by Papadakis and Lawrence (1995) is signif- 
icantly steeper than that of ME curve. On the other hand 
the conclusion of their paper was that power law is not a 
good representation of the power spectrum of LE curve. Vi- 
sual inspection of their Fig. 10 shows that if a single point 
at log{v) ^ —4.2 is removed the mean slope would be again 
consistent with the ME slope and with our result. Perhaps 



there is an additional power at around this frequency as 
well. An extension of the power spectrum towards longer 
timescales based on combined observations on a number of 
satellites operating recently might give the answer. 

4 DISCUSSION 

4.1 Deterministic chaos vs. stochastic signal 

A number of works have boon already devoted to the search 
of low-dimensional chaos in different astronomical objects, 
e.g. Voges, Atmanspacher & Scheingraber (1987), McHardy 
& Czerny (1987), Canizzo & Goodings (1988), Lochner, 
Swank & Szymkowiak (1989), Norris & Matilsky (1989), 
Fiore & Massaro (1989), Cannizzo, Goodings & Mattel 
(1990), Harding, Shinbrot & Cordes (1990), KoUath (1990), 
Krolik, Done & Madejski (1993), LCMH. 

However, the methods used do not provide the unique 
answer to the question unless they are supplemented by 
methods sensitive to the difference between low-dimensional 
deterministic chaos and stochastic phenomenon. One such 
method (nonlinear prediction test of Tsonis & Eisner 1992) 
was applied in this paper whilst two other methods (phase 
randomization method of Theiler 1992 and signal differenti- 
ation method of Provenzale et al. 1992) wore used for analy- 
sis of the optical lightcurve of quasar 3C 345 by Provenzale, 
Vio & Cristiani (1994). 

In both cases the tests indicated stochastic nature of 
the astronomical phenomenon. This shows clearly the im- 
portance of the supplementary tests; any reliable claims of 
deterministic chaos have to pass at least one such test. 

4.2 Consistency of the thermal Comptonization 
model of X-ray emission with the observed 
variability 

On the basis of their results on the variability of NGC 
4051 Papadakis & Lawrence (1995) concluded that thermal 
Comptonization models are ruled out since they always pre- 
dict steeper ME spectra than LE spectra whilst the obser- 
vations give the opposite effect. Since the thermal Comp- 
tonization is currently the most popular mechanism (see e.g. 
Zdziarski et al. 1995) we rediscuss this conclusion within the 
light of our results. 

Since the correlation slope, 2H, calculated for the ME 
and LE curve (see Section 3.4) are the same our conclusion is 
that the character of variability (and the mean slopes of the 
power spectra) are actually the same for both lightcurves. 
In the light of Compton scattering model it implies that the 
multiple travel time within the hot plasma is shorter than 
the timescale for variations of the plcisma parameters. 

As the mean energy photon in LE band is 0.2 kev and 
it is 4 keV in ME band (Papadakis & Lawrence 1995) the 
energy change of the photon by a factor 20 is required to 
change LE photon into ME one. The relative energy change 
is a single scattering is approximately given by 



^ = e(i + e) 



e 



(6) 



where Te is the electron temperature of the thermal plasma. 
The number n of scattering required to interchange LE into 
ME photon can be therefore estimated from the relation 



20. 



(7) 
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This number is about 8 if the temperature of the plasma is 
~ 100 keV and only about 2 if it is ~ 200 keV. 

Therefore, if the plasma temperature is high, the mea- 
sured delay docs not have to bo much smaller than a factor 
few times the light travel time of the emitting region and the 
lack of noticeable delay between the LE and ME lightcurves 
of NGC 4051 (< 60 s, Papadakis & Lawrence 1995) simply 
translates into the size of the emitting region of order of a 
fraction of a Schwarzschild radius adopting the value of the 
mass of the central black hole equal 7 x 10^ Mq. 

It may well be cosistent with the models of relatively 
flat corona, the height being a fraction of the radius and the 
emissivity strongly peaked close to the horizon of the black 
hole. Exceptionally detailed studios of the Ka iron lino in 
Seyfert gala^Ky MCG 6-30-15 (Iwasawa et al. 1996) supports 
similar view, which also strongly suggests that the central 
black hole is rapidly rotating. 

4.3 Constraints for stochastic models from the 
power spectra 

The conclusion from the analysis of the character of X-ray 
light curves of (predominantly) radio-quiet AGN can be ap- 
plied to the detailed models of the X-ray emission. 

The stochastic nature of variability of X-ray radiation 
from active galactic nuclei well agrees with the coronal na- 
ture of emission proposed by Galccv, Rosncr & Vayana 
(1979); De Vries & Kuipers (1992); Haardt, Maraschi & 
Ghisellini(1994); Stern et al. (1995). 

For example, in the solar corona the energy is stored 
in a form of magnetic field loops, and then released during 
reconnections which results in an event of solar flare (sec e.g. 
Demoulin et al. 1993). The efficiency of X-ray production 
is not very high in this case: only some 10% of the flare 
energy is emitted in the form of X-ray radiation (Confield et 
al. 1980). However, the strength of the magnetic field in an 
accretion disk may be considerably higher as the rotational 
velocity is very high (Keplerian) and field amplification may 
operate far more efficiently. 

The stochastic nature of the variability does not proof 
that the magnetic field growing in the result of dynamo ac- 
tion within the disk is the original cause of the variability, as 
some corona models were assumed to be powered by coronal 
accretion itself (e.g. Zycki et al. 1995, Witt et al. 1996) and 
any instabilities growing within such a corona would give 
basically similar effect. 

Unfortunately, most of the models deal mostly with the 
stationary solutions (frequently at a single representative 
radius) and the attention is focused on the radiation spectra 
so at present it is impossible to say whether these models 
are consistent with the variability studies. 

The only model, which can be directly tested against 
the time dependent properties is the clumpy corona model 
developed by Haardt, Maraschi & Ghisellini (1994). 

In the model of Haardt ct al. (1994) the evolution of a 
single clump consists of two phases. The first phase, invis- 
ible and much longer, is the charge time when the energy 
accumulates within a clump due to magnetic field. The sec- 
ond phase is the rapid discharge accompanied with a flash 
of radiation. The charge time in this model depends on the 
radial position of the formed clump {tc ~ r^''^) whilst the 
discharge time is independent from radius. Therefore the 



set of shots consist of the shots with a range of amplitudes 
{Lbiob ~ r^^^) and the same decay times. The power spec- 
tra of such a shot noise are given by a formula (Lehto 1989, 
PapouUs 1984) 

F{u;) = 2n\hld{u;) + ^,_^y^2 - (7) 

where 1/A is the average time interval between shots. Such a 

power spectrum has a shape strongly different from a power 
law with the slope ~ 1.5. More complex coupling between 
the matter and the magnetic field has to be invoked to get 
an agreement with the observed power spectra. 

4.4 EXiture observations 

The method of timeserics analysis presented in this paper 
could be improved by longer data streams (several thou- 
sand points) and a low signal to noise {S/N ~ 1). The 
continguous sections of the data should last at least a time 
~ 10 — 20 times longer than the separation of two significant 
data points. 

5 CONCLUSIONS 

The nonlinear prediction method applied to the EXOSAT 
X-ray lightcurves of Seyfert galaxies indicates that 

(i) there are no signs of the deterministic chaos in the an- 
alyzed X-ray lightcurves of the Seyfert galaxies 

(ii) in four out of eight cases the signal is clearly of stochas- 
tic nature thus generally supporting emission models 
based on multiple active centers 

(iii) in the remaining four cases the variability is not strong 
enough to determine reliably its character but again the 
stochastic interpretation is favored. 

These conclusions apply to the short term variability, 
with the timescales covering from ~ 200 s to ~ 4000 s. 
Long timcscalc phenomena, lasting months - years, may be 
of different nature as the overall spectral changes accom- 
panying large UV brightennings observed in some sources 
are distinct from those observed during the short timescale 
variation (e.g. Done et al. 1995 for NGC 5548; Clavel, 
Wamsteker and Glass 1989 for Fairall 9). 
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APPENDIX A: DATA REQUIREMENTS 

Nonlinear prediction method can be used to study the decay 
of short term correlations in the data and to discriminate 
whether X-ray emission is due to global but highly nonlinear 
behaviour of the emission region or due to multiple flares. 

The presence of noise, however, limits the possibility to 
determine the type of behaviour. The error of prediction 
consists of both the error intrinsic to the dynamical charac- 
ter of the system and the statistical error given as a signal 
to noise ratio S/N. 

The intrinsic error can again be characterized by the 
rate of decay of correlations, (1 — R(T)y^^ and depends on 
the decay law (power law or exponential). The maximum 
time step of prediction is therefore limited by the condition 



[1 - R{T)] + {N/sy = 1. 



(^1) 



It means that the data allow to study the correlations 
up to 



2K 



In 



i (1 - Wsf) 



for exponential decay of correlation and 



{N/sy 



{A3) 



for power law decay of correlations. Here N/S stands for 
noise to signal ratio and A, B, H and K characterize the 
intrinsic source properties. These two formulae allow to es- 
timate the improvement of the data quality with an increase 
of the signal to noise ratio. 
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